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In the context of spinfoam models for quantum gravity, we investigate the asymptotical behavior of 
the {6j}-symbol at next-to-leading order. We compute it analytically and check our results against 
numerical calculations. The {6j}-symbol is the building block of the Ponzano-Regge amplitudes for 
3d quantum gravity, and the present analysis is directly relevant to deriving the quantum corrections 
- - - to gravitational correlations in the spinfoam formalism. 
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J>^i I. INTRODUCTION: SPINFOAMS AND THE {6j}-SYMBOL 

a ■ 

^> . The spinfoam formalism is an attempt to define rigourously a path integral for Quantum Gravity. Spinfoams can 
' be interpreted from several perspectives, as a covariant history formulation for Loop Quantum Gravity describing 
, the evolution of the spin network states, as an improved and quantized version of the Regge calculus for general 
relativity, as a quantization of "almost topological" field theories, as a higher-dimensional generalization of the matrix 
r— i! models generating 2d surfaces. And they have been shown to be related to many other approaches to quantum 
^ ^ gravity. The spinfoam model for 3d quantum gravity is the Ponzano-Regge model which was the first spinfoam 
I ' model ever written. It has been shown to provide a consistent quantization of general relativity in three space-time 
^1 dimensions (for both Riemannian and Lorentzian signatures). The main spinfoam models for 4d quantum gravity are 
the Barrett-Crane model 21 and the more recent EPR-FK-LS family of spinfoam models f3!,'3,'3|. They are related 
to the reformulation of general relativity as a constrained topological BF-theory and were mostly constructed as a 
• discretization of the path integral over space-time geometries. We also point other spinfoam constructions based on 
' another way to discretize the constraints reducing BF-theory to gravity Q, on the quantization of the MacDowell- 
Mansouri action for gravity [7], on the group field theory approach Q and on a generalization of the Barrett-Crane 
model better suited for renormalization [9|. 
. Such spinfoam models provide a description of the quantum geometry of space-time at the Planck scale. The main 
• ' issue is then to extract semi-classical information from the formalism and to show its relation to the more standard 
; perturbative approach to the quantization of general relativity (as a quantum field theory). Solving this question 
' amounts to proving that we recover general relativity in a large scale (or low energy) regime of the spinfoam models 
', and to showing how to compute the quantum corrections to the classical dynamics of the gravitational field. A 
■ proposal to address this problem is the "spinfoam graviton" framework proposed by Rovelli and collaborators p^ . 
It defines the propagators and correlation functions for geometric observables, mainly the area, from which we can 
extract information about the (effective) space-time metric and its (quantum) fluctuations. Most explicit calculations 
$H ' in this framework have been done at the leading order (in the scale parameter) and for the simplest space-time 
triangulation (a single tetrahedron in 3d and a single 4-simplex in 4d). In order to make the link with the standard 
QFT perturbative expansion, we need to be able to push these calculations further and calculate the correlations both 
at higher order ("loop corrections") and for more refined triangulations (smoother boundary state). In the present 
work, we focus on the first aspect: the leading order of the correlations gives the classical propagator of the graviton 
and we would like to compute the higher order (quantum) corrections. Following the lines of [ill. fl^. IT3| . this requires 
understanding the corrections to the asymptotical behavior of the spinfoam vertex amplitude, which is the amplitude 
associated to a single tetrahedron in 3d quantum gravity or to a single 4-simplex in 4d models. This is the basic 
building blocks of spinfoam models, which are then constructed by gluing these spinfoam vertices in some particular 
way in order to describe the whole space-time. In the Ponzano-Regge model, the spinfoam vertex is given by the 
{6j}-symbol from the recoupling theory of the representations of SU(2). The Barrett-Crane model is defined by the 
{10j}-symbol and the more recent models use the EPR or FK vertex amplitudes. The present paper focuses on the 
{6j}-symbol, relevant for 3d quantum gravity. 

They are three basic ways to compute the leading order asymptotics of the {6j}-symbol and show its relation to 
the Regge action for 3d gravity: 
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• Recursion relations fl3]: Using the invariance of the {6j}-symbol under Pachner moves (Biedenharn-Elhott iden- 
tity) or directly its definition as a recouphng coefficient, one can derive a recursion relation for {6j}-symbol. This 
recursion formula is actually very useful for numerical computations, but it can also be approximated at large 
spins by a (second order) differential equation. One then derive the asymptotics from a WKB approximation. 

• Integral formula [isl [iGj: One can write the square of the {6j}-symbol as an integral over four copies of SU(2). 
In the large spin regime, we can use saddle point techniques and one derives the right asymptotics after a careful 
analysis of non-degenerate and degenerate configurations for the saddle points. This is the technique used to 
derive the asymptotics of the Barrett-Crane and EPR-FK vertex amplitudes. 

• Brute-force approximation [17| : One can start from the explicit algebraic formula of the {6j}-symbol as a sum 
over some products of factorials. Using the Stirling formula and after lengthy calculations, we approximate the 
sum by an integral and use saddle point techniques again which lead to the same asymptotics. 

We also point out the more sophisticated and rigorous proof of the asymptotics by Roberts [l^ based on geometric 
quantization, but that also uses an integral formula and saddle point methods. 

The present goal is to push these approaches one step further and derive the first correction to the leading order 
formula. As a first attempt, we focus on the third method and show how to extract the next-to- leading order 
corrections through a brutal approximation of the explicit algebraic expression of the |67'|syinbol. We compare our 
results with the previous calculations for the cases of the equilateral and isosceles tetrahedra [l^ and check the general 
case against numerical calculations. Although the final explicit formula for this next-to-leading order in the general 
case is not particularly pretty, we prove that it is indeed possible to compute it analytically exactly and we show that 
we could extract all orders of the {6j}-symbol using the same procedure. This is a necessary step towards providing 
explicit formulas or procedures to compute all orders of the perturbative expansion (in term of the length scale) of 
the graviton correlations in spinfoam models. 

We can also use the more subtle approach of the recursion relation. This requires a careful analysis of the recursion 
relation and computing the corrections to the WKB approximation [Toj . Or we could use the integral formula 
technique, then one should be particularly careful when dealing with the degenerate contributions to the {6j}-symbol. 



II. THE {6j}-SYMBOL 

The {6j}-symbol is the basic building block of the Ponzano-Regge model which is a state sum model for 3d 
Euclidean gravity formulated as a SU(2) gauge theory. The Ponzano-Regge model is defined over a triangulation of 
space-time: we build the 3d space-time manifold from tetrahedra glued together along their respective triangles and 
edges. We assign an irreducible representation (irreps) of SU(2) to each edge e of the triangulation. These irreps 
are labeled by a half-integer jf, e N/2, the spin, and the dimension of the corresponding representation space is 
given by dj^ = 2je + 1. Each tetrahedron of the triangulation has six edges labeled by six spins je^, ..,^6,3 and we 
associate it with the corresponding {6j}-symbol, which is the unique (non-trivial) SU(2) invariant built from these 
six representations. It is giving by combining four normalized Clebsh-Gordan coefficients corresponding to the four 
triangles of the tetrahedron. Finally, the Ponzano-Regge amplitude for a given colored triangulation is simply given 
by the product of the {6j}-symbols associated to all its tetrahedra. 

Looking more closely at a single tetrahedron, we label its four triangles by / = 0, .., 3. Then each of its six edges is 
labeled by the couple of triangles to which it belongs, (/J) with < / < J < 3. To each edge is attached a SU(2) irrep 

of spin jij, which defines the length of that edge ju + ^ = -y^ (see figure [1]). There are several ways of expressing 




the {6j}-symbol. The basic formula is the Racah's single sum formula which expresses the {6j}-symbol as a sum over 
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some products of factorials (see Appendix Ej) . This is our starting point as in : 

I Joi 302 J03 1 ^ ^A(joi, J02, jo3)A023, J02, J12) A(j23, Ji3, Jo3)AOoi, ii3, J12) V rf u MTT^ / — 7^T 

[ J23 J/13 J12 J ^^^^ n/(i - '"iV- RjiPJ - 

where the vj and pi are given by the following sums: 

VX = 0..3, VK^^jiK, yk^l..3, Pk ^ ^ Uot + jki)- 

I^K i^O,k 

The factors A(joi, jo2, jos) are weights associated to each triangle and are defined by: 

O'oi + jo2 - joaV-Uoi - jo2 + jo3)!(-Joi + jo2 + josY- 



A(joi, J02,i03) 



(joi + J02 + j03 + 1)! 



From this point, in all sums and products throughout this paper, capital indices K will run from to 3 and lower-cases 
indices k will run from 1 to 3. 

We are interested in the large spin expansion of the {6j}-symbol when scaling all the spins homogeneously. Actually 
we will scale the lengths djj j/2 instead of the spins j/j because the structure of the expansion will be simpler (we expect 
an alternation of cosines and sines without any mixing up at all orders as in [T^ 'l and the geometrical interpretation 
(when possible) is expected to be simpler. Then we rescale all djjj by Xdjjj in ([T|), which is equivalent to changing 

jij = ^ - 1 to A%i - i. This gives: 

Xd,,j2 - 1/2 Xd,„j2 - 1/2 Xd,,j2 - 1/2 
Xd,j2 - 1/2 Xd,,j2 - 1/2 Xd,j2 - 1/2 



\/ A(Adjoi , Adjo2 , Xdja^)A{Xdj2s , Xdj^.^ , Xdj^^ ) A(Adj23 , Xdj^^ , Xdjg^)A{Xdjg-^ , Xdj^^ , Xdj^.^ 

A niin p, — 2 . , , 

{t+l)l 



(2) 



..aS-/"'^n.(^-A^. + i)!n,(Ap.-t-2)! 
with the new conventions: 

^K- 2^ Pk = 2^ ^ , 

I^K i^O,k 

\j \j \ _ (2'('^Jcii '^J02 ^ ^ios) ^ 2)' il^'^ioi ^ d.jo2 + ^ioa) ~ 2)' (2'(^^Joi ^ '^Joa + '^joa ) ~ 2)' 
A(Adj„, , Adjo, , Adjo3 j - jj- — ■ 

V 2 ioi ' io2 ^"jos/ 2/' 

The quantity vk gives the perimeter of the triangle K while the p/c's are the perimeters of (non-planar) quadrilaterals. 

III. PERTURBATIVE EXPANSION OF THE 6J-SYMBOL 

In this section, we will give a procedure to obtain the full perturbative expansion of the {6j}-symbol in term of the 
length scale A and we compute explicitly the leading order (Ponzano-Regge formulae) then the next-to-leading order 
analytically. 

A. General procedure 



We give all the necessary formulae to obtain the Ponzano-Regge corrections at any order. But calculations are only 
performed explicitly at the next-to-leading order for a generic {6j}-symbol. We start from equation ([2]). 
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a. First approximation: factorials. The factorial can be expanded in a series: 

/- — /n\« / 1 1 139 571 \ 

nl = V2nn[-) 1^1 + — + -—^- ——^ „ 2488320.4 + ' ' ' j 

In equation there are factorials of the form: nl, (n + 1/2)! and (n — 1/2)!, which are rigourously defined through 
Euler's T function. From ([3]) we deduce asymptotic expansions for (n+1/2)! and (n — 1/2)! (see the details in appendix 
\B\\ . In order to get the next-to- leading order (NLO) in the 1/A expansion of the {6j}-symbol, we replace the factorials 
in equation Q by their respective asymptotic expansion: 

n\ - V2^e("+^) (1 + jL.^ 

(n+i)!^^/2^e("+l)M«)-"(l + _^) (4) 

(n-i)!^V2^e"M«)-n(i __!_). 

Then, equation ^ reads at first order as: 

r Xd,,j2 - 1/2 Xd,„j2 - 1/2 Xd,„j2 - 1/2 1 J_^^Hd,,,) (. _ 
\\d,,j2-l/2 \d,,j2~ll2 \d,,j2-l/2] ' 2A\''^''='' 



0(4) )S. (5) 



.A 

The first factor is given by: 

Kdjij)^^dj,,hd^^^ (6) 

KJ 

with hd — In I ^'^JJJ ~ ^jjif ~t~ djiL){'^ji.j ~^ 'djiK ~ d-jj]:^){djj J — djj^^ + djj^){djj j + djjj^ — djj^) 

2 \ {djjj + djjj^ + djjj^){—djj J + rfj/K + 'djiL){'^ji.j + ^jjjf + d,j j,^){—djj J + djj^^ + djj^) 

where (-fCi) is the opposite side to (/J), that is A' 7^ L and K,L ^ I, J. The second factor is due to the NLO of the 
factorials: 



2Ai 



(7) 



where Aj is the area of triangle /, r^ is the radius of the incircle of triangle / and rj^ is the radius of the excircle to 
the triangle / tangent to the side dj^^^ of the triangle /. Finally, E is a Riemann sum: 



min Pj /2 



2;— max i)//2 



with the prc-factor and the action given by: 

/(x) = iTTX + xln(x) — ^^(2; — vk) — vk) — E](Pj ~ 2;) ln(pj — x). 



1 / n,(P7 - 



13 47 V- 1 10 V- 1 
G{x) = + Tr/ ^ + 13> 



2^X-VK ^Pj-X 



The details of the computation are given in Appendix ICl 

h. Second approximation: Riemann sum. The second approximation consists in replacing the Riemann sum S 
of ^ by an integral. One fc~^ factor of S plays the role of dx. We can then rewrite equation ([5]) as: 

This approximation does not generate any corrections at leading order and at first order. It will nevertheless enter at 
second order in terms in 1/A^. 
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c. Third approximation: saddle point approximation. We have to study an integral of the form / = 
dxg{x)e^-^^^^ where A is a large parameter. The asymptotic expansion of such an integral is given by contri- 
butions around the stationary points of the action / which are points, denoted xo, of the complex plane such that 
f'{xo) = 0. We expand the action f{x) and the function g(x) around the stationary points xq in term of Sx ^ x — xq: 



f{x) = 2^ — ^ — [Sxy = j{xo) + 



3=0 



{5xf+Il\5x) and g{x) = Y1 



— {dxy = g{5x). 



We then expand K{Sx) = g{5x)e^^^o ^^^"^ in power of 5x. Following the standard stationary phase approximation, 
we extend the integration domain to the whole K. The integrals are then "generalized Gaussians" which can easily 
be computed. We group the resulting terms according to their dependence on 1/A, being careful because of the 
function g{x) which depends on 1/A. We recall that g{x) was obtained by replacing the factorials in ([2]) by their series 
expansion and we write g{x) under the general form: 
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Then the complete perturbative expansion of / can be written as: 

"n-l 

A~ 



A/(a;o) 



27r 



-/"(^o)A 



1 + 



(11) 



where 



E[^] 



E 



E[^] (2p-E5 = i'.), ^ i , 

j2 9n-p+, (a^o) -Qr^(a:o) 



(2n)! C,h2 



(a;o)-Q/'^(a;o) 



r2(p+r.)i 



E 



iJ[H2±£)] (2(n+p)-5:5 = i'j) 



(12) 



E 



:ft—p 



11 „ M forp > 1 



The details of the computation are given in Appendix|D] From this expansion and adjusting the first approximation to 
get the proper dependence on A for g and the pre-factors, it is possible to compute analytically the whole asymptotic 
expansion of the {6j'}-symbol. 

Here to get explicitly the next-to-leading order of the {6j'}-symbol asymptotic expansion, we only need the next- 
to-leading order of the 1/A expansion of /, so we cut the previous formulae at n = 1, then 



27r 



V -f'M^ 

with the expansion coefficients given by 



No 



No 



3Ni 



15N2 



-f"M {-f"{xo)? i-nxow 



No^gi{xo), No = 



9o{xo) 



Ni = 



f^'HxMxo) ^ f^^Kxo)go(xo) ^ goM H^^Hxp 



3! 



4! 



V 3! 



We recall that g{x) = e-^(^) (^1- that is: ga{x) = e^(^) and 51(2;) = -^e^("=). Finally, we obtain the 

approximation: 



^^E J -f'%o}^ ^ 



F(xo)+\f(xa) 



[1,1 f Gjxo} _ F"{x„) + {F'{xo)f f(''H^o)+if'^''H^o)F'{xo) _ Hf^''H^o)f \ , n (J-'' 
[^^\ 12 2f"ix„) + 8{f"(xoW 24{f"{xo)r ) ^'-^ KX'', 



(13) 
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This gives us the following expression for the asymptotic expansion of the {6j}-symbol at second order: 

Ad,„,/2 - 1/2 \d,,j2 - 1/2 \d,,j2 - 1/2 
Xd,,j2 - 1/2 \dj,j2 - 1/2 \d^,j2 ~ 1/2 



(14) 



+ A \^ 24 12 2/"(a;o) 8(/"(xo))^ 24(/"(j:o))V ^ ^ ^ P 

where are the stationary points of the phase, i.e. /'(xq) = 0. The next step is to identify these stationary points. 



B. Contributions of the stationary points 

The phase f{x) is an analytical function given by: 



f{x) =mx + x \n{x) - ^[x ~-^]\^[x - ^] - 



VK 

2 



Pj 



K - - - - J 

therefore the stationary points xo satisfy the following equation as shown in [l7| : 

f'{x) = iTT + \n{x) - ^ In (a; - ■5^/2) + ^ In (pj/2 - a;) 

which is equivalent to 

X YliPj -x) = - Yl{x - vk) 

j K 

The previous equation reduces to a quadratic equation Ax"^ — Bx + C = with 

/ \ 

1 



ln(^ 







A = -^pkPi + ^ 



vkvl 



j<l K<L 



I <J,K <L, 
\{I,J)^{K,L) 

~ I 



B = -P1P2P3 + X! ^i^J^K = ^ 

I<J<K 



K 



E 

I <J,K <L, 
J)^{K,L) 



\i<J 



n ^ 



(15) 



(16) 



(17) 



(18) 
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As shown in , the discriminant A = — {B^ — AAC) is given in terms of the djj j by 



A 



16 



E 

/ < J, 

K <L, 
{I,J)^{K,L) 



= 2 



2 
2 



2 



2 



E 

M < N, 
{M,N)^{I,J), 
\{M,N)^{K,L) 



^312 

2 



En^: 



'2 



(19) 





[¥)' 
{¥)' 



1 



(%) 





2 


1 



24(31)2^2 



where V is the volume of the tetrahedron of edge length dj^ jjl. In the following we will focus on the case where 
A > 0, i.e. V^2 > 0, which corresponds to tetrahcdra in flat Euclidean space. The other case A < corresponds to 
tetrahedra admitting an embedding in the 2+ld Minkowski space-time. And so, we get two stationary points: 



2A 



(20) 



The geometrical interpretation of the stationary points is not clear yet. We have shown that A is related to the 
volume of the tetrahedron. B and A are also related to invariant of the tetrahedron: 

B = V — A + 24ycot6i 
^ 2 

where we recall that u/ is the perimeter of the triangle I of the tetrahedron. The angle d is the Brocard angle of the 
tetrahedron. Indeed, iisL^iim. - im.im.ili3. .. iixLi^im. - im^iixLim. are the barycentric coordinates of the 
second Lemoine point of the tetrahedron denoted L. This point is such that the distance from L to the face / of 
the tetrahedron is equal to Ri tan 9 where Ri is the radius of the circumscribed circle of the triangle / and 9 is then 



defined by 



J 



= i2ycot( 



The geometrical significance of the stationary points still has to be understood. However, we can now give the 
explicit form of the leading order and of the next to leading order of the {6j}-symbol. 

Leading order. We first focus on the leading order and on the contribution. This analysis has already been done 
in [l3| and we just recall the main steps and give the notations: 

KJ 



fd =ln 



fd =ln 



The second derivative of / is given by: 



{x^-vo){x+-Vi) 

(x+-Vk)(x+-Vi) 

{Pk-x+)(pi-x+) 



K X-\-—VK 



for ij e {1, • • • ,3} 
for i. A: e {1, ■ • ■ ,3} 

y.^ ^ 

Pj—X+ £C + 



(21) 



= -zVAexp(- ln(x+ Oj-fe - ^+))) 



8 



where we have used the equation ^7} which gives x+ YljiPj — x+) = — JJ^ix-^- — vk)- In the same way, we can 
simphfy F{x+) = — ^ In f a;+ {pj — x+) ) . The exponential piece of /"(a;+) and e^*^^+^ compensate and we get: 



Collecting these different results yields the following contribution of the x+ stationary point: 



-f"i.\,2.X^ exp ) + A/(.+ )) exp (j2^M^h,^^}j 



: exp 



(22) 



\/27rA3\/A 

The same analysis for the x^ contribution yields the same contribution as the previous one with an opposite phase: 



-/ 



.(^^ exp + A/(x_)) exp (j2^h,^^}j 



exp 



(23) 



We must now compute fdj^j which is a complex logarithm. We recall that the principal value of the logarithm is 
defined by Logz :— In \z\ + iArgz. Therefore, we have to compute 3(/dj^ ^ ) = Oij. From (|2ip . we can write that: 



Ooi = Arg(x+ - vq) + Arg(x+ - v^) - ^Arg(pj- - x+) 

Oik = Arg(x+ ~Vk)+ Arg(a;+ - Wj) - Arg(pfe - x+) ~ Arg(pi - x+) 



(24) 



The analysis done in 17] shows that 9jj can be identified as the (exterior) dihedral angles of the tetrahedron. 
Moreover, 



5R(A,J = ln 

A tedious (but interesting) computation shows that: 



(x+-Vk){x+-Vi) 
{Pi-x+){Pk-x+) 



hd- =0 



(25) 



Then, summing the contributions of xj^ and X- we get the leading order of the 6j-symbol: 



Xd,,j2 - 1/2 Xd,„j2 - 1/2 Xd,,j2 - 1/2 \ ko. 
Ad,,3/2-l/2 Ad,,3/2-l/2 Ad,,,/2-l/2 



1 



127rA3y 



cos 



(26) 



(27) 



Xdj 



-9 J J is the Regge action. This is the well-known limit given by Ponzano and Regge [l| and 



where Sr = Y.i<.j — 

which has justified their state sum model for 3d Euclidean gravity where the {6j}-symbol is the spinfoam amplitude 
for a single tetrahedron. 

Next to leading order. The next-to-leading order is then given by the term in in equation (11411 . Using equations 
([7]|B][S]), we rewrite the leading order in terms of x±, vi,pj and A: 



V48^ 



(28) 
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where 



n (^+ - ^^)' + 6n(^+ - E (^+ - ^^)(^+ - *^))^ 

K L^K I K<L 

~H^\{i.pj - - E n (^+ - ^^)' + E4n(^^' - ^+)' 
-(E n (^+ - ^^))(-nfe - ^+)' + E n (^+ - ^^)' - E4n(^' - 

K L^K j K L^K i 

-5(-nfe - + E n(^+ - - E4n(p' - ^^ff\ 

3 K L^K j l^j 



(29) 



Since x± are conjugated to each other, we obviously have = A{x^). Moreover, numerical computations shows 

that 5R(A(a;±, -5/,^^ , A)) = 0, and in particular A{xj^) — —A{x^). This is a priori a non-trivial result to obtain from 
the previous formulas. Nevertheless, we tested it numerically for various choices of spin and it always turned out true. 
Thus we believe that there should be a way to show it analytically. We can then give an explicit formula of the NLO 
of the {6j}-symbol: 



1 



NLO 



: cos 



Sr 



1 



r(A(x+,w/,Pj,A))sin(S'fl + 7r/4). 



(30) 



This result is confirmed by numerical simulations. The plots figure [2] represent numerical simulations of the {6j'}- 
symbol minus its approximation given above pop . Moreover, to enhance the comparison, we have multiplied by A^' ^ 
to see how the coefficient of the next to leading order is approached and we have divided by cos(5'i{ -|-7r/4) (oscillations 
of the next-to-next-to-leading order) to suppress the oscillations; that is we have plotted: 



5nlo = A""/^ 



{6j} - {&j}NLO 

cos(5fl 7r/4) 



(31) 



As expected, the numerical simulations show that this rescaled difference 5mlo goes to as 1/A when A goes to oo. 
Moreover, the data for 5mlo without any oscillation suggest that we correctly divided by cos(S'_r -|- 7r/4) and thus 
the NNLO of the{6j}-symbol should oscillate in cos{Sr + 7r/4). Therefore, this strongly suggest that the asymptotic 
expansion of the {6j}-symbol in term of the length scale A is given by an alternative of cosines and sinus at each order. 
We strongly underline that this is true because we have rescaled the edge lengths dj^ , . If we had instead rescaled the 
spins jij as usually done, we would have found an oscillatory behavior controlled by a mixing of cos and sin at each 
order (as shown explicitly for the case of the isosceles tetrahedron in [l2|)- This suggests that the dj, j are indeed the 
right parameter to consider when studying the semi-classical behavior of the {6j}-symbol. 

The only thing left to do in the present analysis is to provide the NLO coefficient 5(A(x+)) with a geometrical 
interpretation and to show rigourously that 3?(A(a;+)) vanishes. 
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FIG. 2: Plots of the difference Snlo between the {6j}-symbol and its analytical approximation up to NLO. On the left, we 
look at the {6j}-symboI for d\ — 5A, d2 — 7A, da — 9A, d^ — 7A, ds — 9A, de — 9A with the x-coordinate standing for 3A. On the 
right, we've plotted the case di — 15A, d2 = 17A, ds = 19A, di — 19A, ds — 21A, ds = 17A with A running from 60 to 200. 



Finally, we rewrite the approximation up to NLO of the {6j}-symbol in a slightly different manner: 

1 



l + SR + ^^{Aix+)) + 



(32) 
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This shows that the next-to-leading corrections to the {6j}-synibol can be directly considered as corrections to the 
Regge action for (3d) gravity: 

^corrected ^ S R + \^iA{x+)) . 

A 

We point out that an expansion in 1/A with alternating cos and sin could be similarly re-absorbed as corrections 
to the Regge action. This would define in the spinfoam framework the quantum gravity corrections to classical 3d 
gravity due to the fundamental discreteness of the theory. Such correction would enter the gravitational correlations 
(of the "graviton propagator" type) at second order as suggested in [lH . 



IV. SOME PARTICULAR CASES 



A. The equilateral tetrahedron 



For the equilateral tetrahedron, all the edges have the same length: that is V/, J, djjj = d. The tetrahedron with 
edge length d/2 has a volume V — (d/2)'^ \/2/12 and has all equal dihedral angles 9 = arccos(— 1/3). In this case, the 

expressions greatly simplify. For instance, the stationary points are x± = i^^^^^d. Equations (|27p and (|28p reduce 
to: 



NLO 
cqui 



2^/4 / 7r\ 31 

, cos \OR + — ] , 



sin {Sr + ^ 



(33) 



where the Regge action is Sr, — 3d6. The result was already obtained in IT4I. We confirm it by numerical simulations. 
The plot in figl3] gives the equilateral {6j'}-symbol minus its NLO approximation (|33p . Like for the previous plots, we 
have multiplied by A^/^ to see how the coefficient of the next to leading order is approached and we have divided by 
cos(iS'i^ + tt/A) (oscillations of the next to next to leading order) to suppress the oscillations. 



-0.0015 
-0.002 

-0.0025 
-0.003 

-0.0035 ■ 
-0 . 004 



5 ., — &flO 



FIG. 3: Difference between the equilateral {6ji}-symbol and the analytical result (I33p . The x-axis stands for d and d goes from 
200 to 600. 



The isosceles thetrahedron 



We now consider an isosceles tetrahedron that is a tetrahedron which has two opposite edges of length equal to ^ 
and ^ and the remaining four edges of the same length equal to | (see figure S]). The volume of the tetrahedron is: 



28(3!)- 



;dldl (4d2 -dj- dl) 



and the dihedral angles are: 



9 — arccos 



—did2 



9i 2 = 2 arccos 



d2S 



Ud^ - dl 2 



Once again, equations ([27)1 and (|28|) simplify and we get 
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FIG. 4: The isosceles tetrahedron 



1 



\/l27rT/A3 



cos 



^ ^ 7r\ F{d, di,d2) 
" 4/ 2WXVl2nV)^ 



Sn + - 



(34) 



where F{d, di, c?2) 



76Sd^{d^-dl-dl)+736d*dldl+24:0d'^{d'l+d^)-176d'^dldl{dl+dl)-24:d'^{dl+dl)+Wdldl{df+^ , 

^6*2- Let us point out that the volume increases as A'^ while F goes as A^, so that 



Regge action Sr = 2de + ^Bi 
the NLO scales properly as A~^/^. 

This reproduces the result previously obtained in Jv^. We can easily check that this reduces to the previous 
equilateral case when di = d2 = d and we further confirm it by numerical simulations. The plots in figure [5] represents 
numerical simulations of an isosceles {6j}-symbol minus the analytical formula (|34p . Like for the previous plot, we 
have multiplied the data by A^/^ to see how the coefficient of the NNLO order is approached and we have divided by 
cos{Sii + 7r/4) (NNLO oscillations) to suppress the oscillations. Finally, the geometrical interpretation of the term 
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FIG. 5: Differences between isosceles {6j}-symbofs and their anafyticaf approximation (|34|) . The x-axis stands for d with A 
goes from 200 to 600. On the left hand side, we consider isosceles tetrahedra with di = 3A, d2 = 3A, d = 7A. On the right hand 
side, we've plotted the case d\ = 9A, d^ = 3A, d — 21A. 



F[d^d\,d-2) remains to be understood. If we can't provide it with a geometrical meaning, there is little hope to 
interpret the NLO coefficient 9(A(x+)) in the generic case. Nevertheless, we give a more compact expression for the 
denominator of F: 

96 (4^2 - dX) (4^2 - £^ (4^2 ^d\- dl) = 96^—^. (35) 

We still need to express the numerator of F in term of geometrical objects. For instance, we could express it in term 
of d^, (4(i^ — di){4:d'^ ~ c^i) {4d^ — d\ — d^)^ which would provide a formula in term of the volume and the dihedral 
angles. Nevertheless, we haven't been able to find such a useful rewriting of this NLO coefficient. 



Conclusion 



We investigated the asymptotical behavior of the {6j}-symbol. Starting from its expression as a (finite) sum over 
(half-)integers of algebraic combinations of factorials, we followed the footsteps of [l3| and showed that one can derive 
systematically the corrections to the leading order formula at any order. The method relies on three steps. First, 
we use the Stirling formula (with the appropriate corrections) to approximate the factorials. Second, we consider the 
sum as a Riemann sum and approximate it by an integral (over the real line). Finally, we perform a saddle point 
approximation to compute the behavior of the {6j}-symbol for (homogeneously) large spins. 

Using this framework, we showed that we recover an oscillating leading order (LO) with frequency given by the 
Regge action as is well-known and was already proved in jl7i] . Then we computed analytically the next-to- leading 
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(NLO) corrections. The formula that we obtain is exphcit, although not compact, and we could not interpret it 
geometrically in a clear way. Nevertheless, we performed two simple checks. First, we checked that our complicated 
formula reduces to the known expression for the NLO for isosceles tetrahedra [l2| . Second, we checked it numerically 
in various cases and found a perfect fit. These numerical simulations also confirmed that the NLO is a -j-phase shift 
with respect to the LO (the NLO is given by a sin instead of a cos) and that the NNLO is back in phase with the LO 
(again a cos), which confirms our expectation of an alternating asymptotical series in cos +i sin +p- cos sin + . . . . 

We point out that we computed in details the corrections due to the Stirling formula and to the saddle point 
approximation. However we didn't study the Riemann sum approximation. It does not contribute to the LO and 
NLO. It will only enter at the level of the NNLO. 

This work is mainly technical and can be applied to the computation of gravitational correlations for 3d quantum 
gravity following [T^ . [ill . [T^ . [T3| . It will enter the quantum corrections to the propagator/correlations at second order, 
as was shown in 11]. Indeed, the first order corrections are derived from the path integral of the Regge action, while 
the deviations from the Regge action as computed here enter at second order (as two- loop corrections). From this 
perspective, this NLO of the {6j}-symbol describes the leading order deviation of quantum gravity with respect to 
the classical gravity. 

Beyond the technicality of the paper, our purpose was to show that computing such corrections is indeed possible 
(although it does lead to complicated expressions) and that similar methods could be used for 4d spinfoam gravity. 
Although these methods allow straightforward (but lengthy) analytical calculations, which might be handled by a 
computer program, their drawback is the loss of the geometrical meaning of the expressions obtained. An alternative 
way to proceed is to use the exact recursion relations satisfied by the {6j}-symbol (see 14]) and other spinfoam 
amplitudes (see [l^l) to probe the asymptotic behavior and the induced corrections of the correlations. This is left to 
future investigation p^ . 
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APPENDIX A: THE {6j}-SYMBOL - RECOUPLING THEORY 

The {6j}-symbol is a real number and it is obtained by combining four normalized Clebsh-Gordan coefficients along 
the six edges of a tetrahedron, with edge lengths given by jjj + | = (0 < / < J < 3). We usually express the 




FIG. 6; Tetrahedron: edge lengths are given by hj = 



6j-symbol in term of the Wigner 3j-symbols : 

- / ^y'oi+ioa+joi— Qoi— Q03 — aoi ( -^"l -^12 Jl3 \ 

" \ aoi ai2 -ai3 J 



(Al) 

Jl3 J23 Jos \ / J03 J02 JOl \ / 302 J23 Jl2 



ai3 Oi23 "03 / \ ao3 ao2 -OiQi J \ ao2 0123 ai2 
The Wigner 3j symbols are very simply related to the Clebsh-Gordan coefficients < joiji2Q!oiQ^i2bi3Q^i3 > by: 

< JoiJi2aoiai2|ji3ai3 >= (-l)^--^-+"-(2ji3 + 1/2)'/^ ( ^"^ ) , 
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And Racah gave a general formulae for the Clebsh-Gordan coefficient: 
< ioiil2Q!oiQ!l2|il3Q!l3 > = S{aoi + ai2, Q;i3)A(joiii2ii3) 

\/(2ii3 + l)(ioi +aoi)!(joi -Q;oi)!(ii2 +ai2)!(ii2 -Q;i2)!(ii3 +ai3)!(ii3 - ais)! 

y- (-1)" 

where A(.7oi.,7i2,.n3) = From these, Racah gave a tensorial formulae for the 

6j-symbol, the Racah's single sum formulae: 



JOl J02 J03 
^23 jl3 jl2 



= V A(ioi , io2, ios) A(j23, io2, ii2)A(j23, ii3, ios) A(joi, il3, il2 

nin pj 

E(-i) 



minp, ^ (A2) 



max vj 



Uiit-viy-UM-ty- 



with vk = '^jiK Vii" e {0, • • • ,3} and = ^ (joi + jw) V/s e {1 • • • 3}. 

APPENDIX B: FACTORIALS 

The factorial n! is defined for a positive integer n as: 

n! = n(n-l)---2-l=r(n+l), 

where r(n) is the gamma function for integers n. This definition is generalized to non- integer values. Using the 
identities for the F function, we write explicitly the values for half-integers: 

/ l^■ /— /In, / In. \/7r(2n)! , 1,, -v/tt A/7r(2n + l)! 

where n!! is the double factorial : 

r n • (n - 2) • • • 5 • 3 • 1 if n > odd, 
n!! = < n • (n - 2) • • • 6 • 4 • 2 if n > even, 
(1 if n = -1 or 0. 

Using the asymptotic expansion of a large factorial n! ~ \/27rn (f )" (l + 15;^ = asS;^ ~ sialon-^ ~ 2488320n^ " ' )' 
can get an asymptotic expansion for: 



in + 1/2)1 ^ V2^e("+i) (1 + _L) (1 + 



139 571 



2(2n) ^ 288(2n)2 51840(2n)3 2488320(2n) 

/i 1 1 I 139 ' I 571 _ \ 

\^ 12n 288n^ 51840n^ 2488320n* "')' 



(Bl) 



fr, — iV ^ ,/9^^"ln(")-" ( 1 J- 11 I 1 139 571 , 

ya ^j. \ Z.TIC I ^ 12(2n) 288(2n)^ 51840(2n);^ 2488320(2n)* ' 

( \ 1 1 I 139 I 571 _ ^ 

12ri 288n^ 51840n^ 2488320r)J ' ' ' / ' 



or more simply, at the next-to-leading order: 

n!~x/2^(^r(l + ^), 

(n + i)! ~ ^e^n+l)\n{n)-n (l + ^) , (B2) 
(n-i)!~ V2^e"i"(")-"(l-2l^). 
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APPENDIX C: FIRST APPROXIMATION: FACTORIALS — > NEXT TO LEADING ORDER OF THE 

STIRLING FORMULA 

In this section, all computations are done at the next-to-leading order. We replace the factorials in equation ^ by 
their respective asymptotic expansion. 

• Then, a typical triangle coefhcient: 

A/\ , \, , \ _ ('g^'^Joi ~^ ^302 ~ '^303) ~ 2)' il^'^joi ~ '^jo2 + ^joa) ~ 2)" il^^'^Joi ^ ^J02 + ^jos) ~ 2)' 

^[^"■joi 7 •^"i02 : ^"303 ) — 7X71 7~3 7~1 N IN, 

V 2 v%oi + '^302 + "io3 J 2 / • 

will be 

AfAd,- Ad, Xd, ) = 27r[e"^(''joi+'^™2+'^jo3) ln[|(f'joi+'^jo2+'im3)l + |(dmi+<^io2+'^TO3) f l n 1 , 

("^301 +^j02 ~''j03 ) '"[■5 ("^301 +''j02 ~''j03 ('^301 +''j02 "'^303 ) | 1 — . i- ^ | 

g"2 (^301 ~'^302 ^^303 ) ^^[^ (*^301 ~*^302 ^"^303 )] ~ "2 (^301 ^^302 ^^303 ^ ( 1 



g-j ( — djjjj +djQ2 +'^303 ) '"It (~'^301 +"^302 +'^303 )1 ~ T (""^301 +"^302 +'^303 ) ( 1 , 3— ■ ^ | 

which simplifies 

A{Xdj,-,^, Xdj„2, Xdjg^) = 27re^l'^"'^J01+'^3n2+'*303)l"(-''301+'*302+<^303) + ('i301-'i302+'i303 )ln('i301-''302 +'^303)] 

g- ! [(''301 +'^302 -^^303 ) 1"('*301 +''302 -''303 ) + ('*301 +^^302 +'^303 ) Mdjol +'^302 +'^303 )1 (CI) 
[1 _ J-( 1 + 1 + 1 1 )1 

L 12A ^ -dj„j+djj,2+(ij„3 <ij(,i-d-jo2+d3(,3 dj^-^+d^g^-d-jaa djg-^+dj^^+djg^ n- 

The factor a/ A{Xdj„, , Adj^a , >^dj„^)A{Xdj2.^ , Ad^o^ , Adj^ J A(Adj23 , Xdj^^ , Xdj„^)A{Xdj„, , Ad^.g , Adj^^ ) in equation 
([2]) can then easily be put into the form: 

(2^)2e4'^('^3..)A__l_iy(d^.^^)') (C2) 



KJ 

with /irf = i In '■'^^ij^'^^iK+dj,L)(djjj+dj^^-dj^^)(djjj~d,j^+djj^){djjj+djj^~d-jjL) 



where 



■3/j 2 \{d-jij+djj^+djj^){-djjj+djjj^+djj^){djjj+d.i 

K ^ L and K,L ^ I,. J ^ ' 

H{dj,,) = 2^-^^-: '^here X e {0, - • • ,3} and j G {!,• • • ,3} 

and we recall that {i^ = ^ ^ Vif e {0, • • • , 3}, = ^ (^J"- Vfc e {1 • • • 3}. 

We now replace the factorials in the sum of ^ by their approximations and we change of variables: t — Ax: 
, _ , (Ax + l)(Ax)!n,(Afe-x))n,(Afe-x)-l) 



n/(A(a; - vi) + 3/2)(A(x ~ 5,) -f 1/2)! n,(Afe - x))! 

mill pj 

a;— max vi 
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where 

Gi{x) = inXx + 31nA + Ina:: + 2^\n{pj - x) - ^ ln(a: - vi) + {Xx + l/2)(ln.T + InA) - Aa: + ^A(.t - vi) 

-^(A(x - w/) + l)(ln A + ln(a; - vi)) - ^(A(pj ~ x) + l/2)(lnA + ln(pj - x)) + ^Afe- - x) 
I j j 

(C5) 

which can be simpHfied using the fact that vi — Pj- 

V + xhix ~ ^^(a; — vi) \n{x — vj) — ^^(Pj — x) ln(pj — x) 



1 x^YlAPi-x)^ 
Gi (x) = -2 In A + - In ''p"^' , ; + A 

2 lljix - vi)'^ 



and 



G2{x) 



l-l^-^+y 47 ^ 13 /I 



(C6) 



(C7) 



Moreover, 



where 



/(x) = iTTx + X \n{x) — ^(a; — vk) ln(a; — w/) — ^(Pj — a;) 



F(a;) = ihJ^iI>^ 



2 1^ 0/(^-1'/) 

Then the sum can be approximated by: 



(C9) 



mm Pj . / \ \ 

^<"-)=(sk^ i: .»«->™ l-yi^GW+O i (CIO) 



a;— max vi 

where 



G{x) = -H + + E (Cll) 

X 24(a; — vk) Pj — x 



APPENDIX D: THIRD APPROXIMATION: THE STATIONARY PHASE METHOD 

We are interested in the 1/A expansion of the integral: 



■/max 57/2 V -L^^ V-^ 



We do not give here the proof of the whole expansion (equation (fTTjl ) because of the heavy formalism but we directly 
prove the next to leading order formula (equation p^ ): the procedure is the same but the computations are easier. The 
asymptotic expansion of such an integral is given by contributions around the stationary points of the phase denoted xq. 
We expand the phase f{x) around the stationary points xq at fourth order and the hmction g{x) = e^'^^'^ (l — y^G(x)) 
at second order and we extend the integration to infinity. 



+00 



I^y d{5x) ig{xo) + g'{xn)5x + y{xa)(5xf) eKf^^-^+'^fMi^^?) 

^J-^ ' (Dl) 

\ + A (i/(3)(xo)(fe)3 + yi^){x,){5xf) + f {yi^){x,){5xff + 0{\^)] 
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where in our case, g{x) = e^^^'' (l — j^Gix)) and then the integration are "generahzed" Gaussians: 



Xo 



/•+00 „ 



r-\-OC -f"( 1 



+A 



(£o) , ri(£n) p/ 



2 V 3! 



2 



which can easily be computed: 



2zr pi=^(a:o)+A/(2:o) 



/"(2:o)A' 



L 1 / Gixg) F"{xo) + {F'{xo)f , /W(xo)+4/(^)(a;o)f'(xo) SC/^fW)" A , n M 
+ A 12 27^) + 8(/"(xo))^ 24{f"(xo}r )+^\\) 



(D2) 



(D3) 
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